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ABSTRACT 



Unambiguous accurate direction finding results are 
obtained by the use of large stationary antenna systems 
at the lower frequencies and by the use of rotatable antennas 
at the higher frequencies. While these techniques result in 
good data they suffer from the fact that accuracy is general- 
ly a direct function of the antenna size, that is, accuracy 
increases as the antenna size increases. 

A system is proposed which encodes the outputs of 
several simple omnidirectional antennas and processes this 
modified signal to provide angle-of-arrival data. The 
processor not only provides the angle-of-arrival data but 
the received signal is preserved as a simultaneous output 
with the bearing data. 

The coding is accomplished utilizing complementary 
sequences and a derivative of these sequences called orthog- 
onal complementary sequences. The technique is described 
in detail and a digital computer model was developed to 
support the theory. 
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I. INTRODUCTION 



The f undainental method for determining the location of 
a radio frequency (RF) emitter lies in defining a set of 
independent curves in the required two or three dimensional 
space with each curve containing the emitter location as 
one of the sets of points of the curves . The intersection 
of these curves will locate the emitter . 

The function of a direction finding (DF) system is to 
determine the unambiguous direction-of-arrival (DOA) of an 
incoming signal relative to a fixed set of coordinates at 
the DF site. This DF system would in practice be one of 
several such systems in a coordinated network comprising 
an emitter locating system. Although absolute accuracy of 
DOA information is paramount the ability to determine this 
information in real, or near real, time is generally 
critical . 

Therefore in real applications the ability to determine 
the angle-of-arrival (AOA) of a signal in a timely manner 
is the function of a DF system which is utilized in a 
coordinated search for emitter locations. 

A. CONVENTIONAL DIRECTION FINDING TECHNIQUES 

There are four basic techniques for generating the 
data required to determine the AOA, they are: 

1) Relative amplitude. 

2) Relative phase. 

3) Time-of-arrival (TOA) . 

4) Doppler. 
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1 • Relative Amplitude Techniques 

These techniques utilize directional antennas and 
are quite common below 200-300 MHz. They can be catego- 
rized into three general methods. 

a. Rotatable Antenna Systems 

In this system the entire antenna and therefore 
the associated directional pattern is rotated about a local 
axis. The line of bearing to the emitter is determined by 
observing the angular position of the antenna when the 
response is either at a maximum or at a minimum value. 

b. Rotatable Pattern-Fixed Antenna Systems 

Systems in which the antenna is too large to 

be physically rotated at the desired rate operate by 
rotating the associated directional pattern of the antenna 
or array. This rotation is usually accomplished electro- 
mechanically . The line of bearing is determined as in the 
rotatable antenna system. 

c. Instantaneous Systems 

This system employs two or more antennas and 
receivers in a direct amplitude comparison of the signal 
from the overlapping directional patterns of the antennas. 
The line of bearing is determined by comparison of the 
relative signal strength in each directional pattern. 

2 , Relative Phase Techniques 

These systems utilize non-directional antennas to 
obtain directional AOA data. This technique relies on the 
fact that a wavefront incident upon an antenna system will 



9 



be received at the individual antennas in the system with 
differing phase components. The comparison of the phases 
will be a direct indication of the AOA of the incoming 
wavefront (signal) . 

3 . Time of Arrival 

This technique requires the time measurement of a 
given modulation event at two or more antenna locations. 

The time dif ference (s) of the event at the antennas is/are 
then utilized to determine the AOA. This technique is 
similar to the relative phase technique except that the 
parameter measured is time vice phase. The requirement 
for a well defined modulation event is of importance in 
this method. 

4 , Doppler 

The doppler technique relies on the principle that 
a receiving system in motion relative to the incident wave- 
front of a signal will produce an output in the receiving 
system which is a modulated representation of the original 
signal. The extraction of this frequency modulation 
(doppler shift) information is then processed to determine 
the AOA data [l] . 

B. DISADVANTAGES OF CONVENTIONAL DF SYSTEMS 

All of the above techniques have certain disadvantages. 
The relative amplitude method produces accuracy in direct 
relation to the directional ability of the antenna (s). 

That is a large many wavelength antenna system will produce 
a narrow beam or beams and therefore the resolving ability 
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of the system will increase resulting in an overall increase 
in bearing accuracy. The major drawback to this technique 
lies in the ability to produce the required antenna within 
the physical limitations governed by the application. 
Additionally if the system must be physically rotated then 
there is an additional design burden imposed on the ability 
to obtain the required AOA accuracy. 

Relative phase techniques lead to ambiguous data unless 
the separation between the two antenna is small; less than 
half a wavelength. Although this spacing eliminates ambig- 
uous data due to separation it is in opposition to the 
requirement for accuracy; accuracy requires wide antenna 
separation. 

Time of arrival methods, while quite good above 1000 
MHz especially for pulse signals, rapidly deteriorates 
below approximately 300 MHz due to difficulties in locating 
the same modulation event precisely at each antenna 
location. 

Doppler techniques require a circular array of antenna 
elements with a specified diameter of several wavelengths. 
The physical size of this antenna system often precludes 
its use in many applications. 

C. A NEW DIRECTION FINDING TECHNIQUE 

An approach to providing the required direction finding 
data utilizing a phase coding scheme similar to that of 
intrapulse modulation is proposed. 
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This phase coding technique is similar in philosophy 
to that attributed to P. M. Woodward for active radar 
applications [2] and described by G. J. Sieren in his 
Master’s thesis [3]. 

Basically this technique consists of coding the outputs 
of several antennas, summing the outputs and then simulta- 
neously producing an pseudo autocorrelation and a cross- 
correlation output. The antenna outputs are coded 
utilizing a complementary sequence. The pseudo auto- 
correlation hereafter called autocorrelation, is accomplished 
in the usual manner and the cross-correlation is accomplished 
utilizing a totally orthogonal complementary sequence 
generated from the original complementary sequence. The 
output of the matched cross-correlator will be zero and 
will be unique for a specified angle-of-arrival of the 
signal. Additionally the output of the matched autocor- 
relator will be replica of the incoming signal but of an 
amplitude equal to twice the code length. The angular 
resolution will only be limited by the ability to physically 
produce differing sets of delays in the correlators. 

The next chapter will discuss phase coding in general 
providing a brief description of some codes and then a 
description of complementary sequences (codes) including 
totally orthogonal complementary sequences (codes) will 
be given. 

With this background the next chapter will describe 
the proposed direction finding system in detail including 
results obtained from a digital computer simulation. 
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The last chapter contains a summary, conclusions and 
suggestions for further investigation. 

The appendix describes the digital computer model 
including definition of terms, uses and the complete 
computer source program. The computer program was written 
in FORTRAN IV, Level G and was processed on an IBM 360/67. 
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II. PHASE CODING 



Phase coding techniques and various phase codes will 
now be described in order to provide a basis for the fol- 
lowing chapter which describes the direction finding 
technique utilizing phase coding, 

A signal can be coded (divided) into n equal segments 
of duration T. If the only allowable phase states in each 
segment relative to a reference are zero degrees and one 
hundred eighty degrees then the coding is termed binary 
phase coding. If m possible phases are allowed in any one 
segment then the coding is called polyphase coding. Phase 
codes can be further subdivided into two general classes 
depending on the length and the periodicity of the code. 

The codes can be continuous or confined to a finite time 
segment. Continuous codes find wide use in communication 
systems and in CW radars. Coding confined to a finite time 
segment are generally employed in pulse radar applications. 
If the code is repetitive it is termed periodic and if not 
it is termed aperiodic coding. A continuous transmission 
which has been random binary coded is a special case of 
aperiodic coding. Coded segments, coded pulses and coded 
words are in general synonymous terms. 

A. POLYPHASE CODES 

A random binary coding can be described as dividing 
the carrier into equal segments and assigning a plus (+) or 
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a minus (-) value to each segment in a random sequence. 

The plus and minus values refer to no phase reversal (+) 
and to 180 degrees phase reversal (-) . This representa- 

tion is shown in figure 2-la. The amplitude of the kth 
element in figure 2-la is aj^. The autocorrelation function 
of this waveform can be expressed [4] in discrete form as: 

n 

$(m) = Z *-(n-l) f. m ^ (n-1) 

k=l 

The symbol m is an integer index which steps over the domain 
shown. Each step corresponds to a shift in of one 

segment of duration T, Although this description of the 
autocorrelation function is not quite accurate this discrete 
procedure simplifies the calculation and sketching of the 
autocorrelation function. The autocorrelation function is 
always maximum at t (tau) equal zero and is equal to n. 

For the condition of randomly chosen polarities of a^. and 
for n large then the autocorrelation function near the 
origin is triangular as shown in figure 2-lb. The power 
spectrum of a random sequence, according to the Weiner- 
Khintchine theorem, is the Fourier transform of its 
autocorrelation function [5] and therefore the shape of 
the power spectrum is approximately that for a single 
segment of length T. 

The value of the autocorrelation function away from 
the origin on the time axis can be shown by utilizing the 
discrete representation of the autocorrelation function 
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described above. Assuming a code as follows, +++--+ 
and an offset m equal to one segment, then we have: 

+ + + — — + — 

" +^-+ + - - + - 

1 1-1 1 - 1-1 

Therefore the value of the autocorrelation function equals 
zero which is the sum of the products. In forming the 
products it is assumed that like signs produce a +1 and 
unlike signs a -1. If m=2 then the application of the 

above produces a value of $ (m=2) = -1. 

These sidelobes, the value of the autocorrelation 
function away from the origin are sometimes known as range 
sidelobes. The reduction of the height of these range 
sidelobes is of paramount importance in radar applications 
The coding of the pulse in a radar system results in an 
increase in the range resolution due to the pseudo short 
(width=T) pulse while still retaining the high average 
power advantage of the long pulse (width=nT) . For a long 
random code the value of the range sidelobes will vary 
between ± 0.77n with a probability of 0.5 [4]. 

Polyphase codes can be generated for the periodic case 
which will have zero periodic correlation except for a 
shift equal to zero [6]. 

The correlation function will have the familiar 
triangular shape for a rectangular pulse input. A modula- 
tion code of length n can be expressed as a sequence 
{aQ,ai ,a2 , . . . f^n-l^ • periodic correlation function 

is defined to be the sequence {xq ,X j^ ,X 2 , . . . } . The 
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general function is as shown below. 
n-1 

= E ^k+l^k* Where: * denotes the 

k=0 complex conjugate. 

The definition of the autocorrelation function will result 
in values for integer shifts. The sequence { aQ , . . . , aj^_ 3 ^} 

o 

can be generated for any given code length n=p"^. Where p 
is a prime number greater than unity. The operations are 
as follows. 

1. For a given code length n pick the prime number p. 

2. Find all the pth roots of unity. That is roots of 

the form e“J for ^ ^ ^ P“l» 

3. Form the sequences of length p as follows: 

1 , 1 , ... 1 

R^, _ R^P-1 

2 P“1 

^p-l'^p-1' ••• ^p-i 

The p roots of unity are: 1, Ri , R« , R t- 

X z P”1 

4. The code sequence is formed by writing in serial order 
column one, column two, column three, etc. 

This code thus formed has a correlation function Xj^ which 

o 2 

has a zero value everywhere except for i=0, , 2p , 3p^, 

etc. Although this is a good code (all zeros except main 
peaks) it suffers from the requirement that there be p 
different phase shifts in the code. 



1, 

1, 



1, 
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The restriction that p be a prime number has been 
proven [7] by R. I. Frank and S. A. Zadoff and confirmed 
by R. C. Heimiller not to be a prerequisite. Heimiller 
had originally specified p to be prime. 

A patent by S. A. Zadoff and W. Abourezk [8] utilized 
a basic form of the Frank code sequence to describe a matrix 
as follows: 



1, 2, 3, 4, 

2, 4, 6, 8, 

3, 6, 9, 12, 

4, 8, 12, 16, 



P 

2p 

3p 

4p 



. . . p 



2 



There are p rows and p columns. The numerals represent 
the multiplying coefficients of a basic phase angle. The 
basic phase angle is defined as 2iTr/p where r is relatively 
prime to p and r and p are integers. Letting R equal the 
primitive pth root of unity then the following matrix is 
formed. 



R r2 r3 ... rP 

R^ R^ r6 ... r2p 

R^ R^ R^ ... R^P 



rP 



2 
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Where the exponents of R are given by the first matrix. 
This is just a permutation of Heimiller's derivation. 
Heimiller's proof [6] of the good properties of the auto- 
correlation function applies equally as well to the above 
for p not prime. A special case is p=2 the binary case. 
For p=2 the 2x2 matrix of the pxp matrix is used as 
follows : 



2x2 matrix 



then 



and 



1 2 
2 4 



R r2 
r2 r4 



where r = e“3 <2itr/p) 
= e-j(2Trr/2) 
= e“3 (irr) 



e-jir 


e-j2ir 




-1 


1 


- j 2 tt 
e 


-i 4 tt 
e 




1 


1 



Therefore the code sequence is: -1, 1, 1, 1. 

Polyphase codes with good non-periodic properties can 
be generated by a permutation of the previously given 
matrix [9]. This matrix is: 



0 0 0 0 

0 12 3 

0 2 4 6 

0 3 6 9 

6 



0 

(p-1) 

2(p-l) 

3 (p-1) 

• 2 

(p-1) 
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The phase angle is still 2-nr/p with r relatively prime to 
p. For clarity the periodic philosophy will be restated. 

1. The code sequence is formed by writing the matrix a 
column at a time. 

2. The entire code is periodic that is the matrix is 
repeated. 

3. The autocorrelation function is everywhere zero 
except at multiples of the code period where it is 
a single maximum. 

A non-periodic code can be formed by selecting a starting 
point, running through the matrix once and stopping. The 
assumption in doing this is that there is zero amplitude 
before and after the sequence. There will be more than p2 
starting points. The position in the matrix where the 
phase changes are the least drastic will be the starting 
point. This position seen by inspection is the upper 
left hand corner of the above matrix. 

For a code of length n=16 the basic phase angle is 
euqal to 2iTr/p where r=l is the primitive pth root and 
therefore the general element in the matrix is e”^ (2ir/p) ^ 
The non-periodic exponent matrix is : 

0 0 0 0 

0 12 3 

0 2 4 6 

0 3 6 9 

and the non-periodic phase matrix is : 
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1111 
1 “1 



1-1 1-1 
1 e-j^/2 



The code sequence is obtained by writing in serial form 
column one, column two, column three and column four as 
follows : 



This code has four phase states: zero, 90, 180 and 

270 degrees. The autocorrelation function for n=16 shows 
[ 9 ] that the main peak is equal to 16 and the maximum 
sidelobe is 1.4. This value for the maximum sidelobe is 
found by the vector addition of p/2 unit vectors at a 
phase angle of 2ir/p for p=even or the vector sum of 
(p+l)/2 unit vectors at the phase angle 27 r/p for p=odd. 

There is a class of binary coded pulses whose auto- 
correlation over one period (finite time autocorrelation) 
has only two non-zero sidelobe levels. These codes are 
called perfect words and include Barker codes. The known 
Barker codes are the following. 




- 1 , 



Code Length 



Code 



2 



+ 



3 



+ + 



4 



+ + + - and + + - + 
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Code Length 



Code 



5 + + + - + 

7 + + + -- + - 

11 
13 

The autocorrelation function for these codes is: 

fn if m = 0 

<!> (m) = < ' 

(zero or ± 1 if m 0 

B. COMPLEMENTARY SEQUENCES 

All these polyphase codes including the binary phase 
codes have the property that the autocorrelation function 
is maximum at t (tau) equal to zero but all including the 
perfect words have autocorrelation sidelobes which have 
non-zero values. 

There is a series of binary elements known as comple- 
mentary series which eliminates this non-zero sidelobe 
problem. 

A complementary series is defined by M. J. E. Golay 
[ 10 ] to be a pair of equally long finite sequences of two 
kinds of elements which have the property that the number 
of pairs of like elements with any one given separation in 
one series is equal to the number of pairs of unlike 
elements in the other series with the same separation. If 
the two kinds of elements in the series are +1 and -1 then 
the sum of their autocorrelation functions is everywhere 
zero except for the center term. 
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Assume a complementary sequence pair of length equal 
to four as follows with the two elements being +1 and -1. 



Code A +1, +1, +1, -1 

Code B +1/ +1/ “1/ +1 



Code A has exactly two pairs of adjacent like elements 
(+1, +1, ... +1, +1) while code B has exactly two pairs 

of adjacent unlike elements (+1, -1, ... -1, +1). Addi- 
tionally for elements spaced by one- intervening element 
code A has exactly one pair of like elements (+1, +1) and 
code B has exactly one pair of unlike elements (+1, -1) . 

These like and unlike combinations continue for all N-1 
possible separations where N equals the number of elements 
in the series. 

Golay has shown [lO] that N must always be an even 
number and must be the sum of two squares. Therefore the 
possible values for a code length up to for example 50 
are 2, 4 , 8 , 10, 16, 18, 20, 26, 34, 36, 40, and 50. 

Kernals are basic code lengths which cannot be de- 
composed into shorter lengths. Kernals exist for codes 
of length 2, 10 and 26 [l0,ll], Golay has shown [lO] that 
there are 2^ pairs of complementary series (some identical) 
which can be generated by the following six operations on 
a given pair of a complementary series. 

1. Interchanging the series. 

2. Reversing the first series. 

3. Reversing the second series. 

4. Altering the first series. 

5. Altering the second series. 

6. Altering the elements of the even order of each series. 
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Altering a series embodies interchanging the kinds of 
elements. For example the element +1 becomes -1 and 
conversely the element -1 becomes +1. 

As an example of the above operation consider the 
following example. The original code is: +1, +1, +1, -1 

and +1, +1, -1, +1 and the operations are the following. 



1. Interchange +1, +1, -1, +1 and +1, +1, +1, -1 

2. Reverse-first -1, +1, +1, +1 and +1, +1, -1, +1 

3. Reverse-second .... +1, +1, +1, -1 and +1, -1, +1, +1 

4. Alter=first -1, -1, -1, +1 and +1, +1, -1, +1 

5. Alter-second +1, +1, +1, -1 and -1, -1, +1, -1 



6. Alter-even order .. +1, -1, +1, +1 and +1, -1, -1, -1 
All of the above codes and the remainder of the sixty- 
four combinations will maintain the properties of comple- 
mentary series. 

Golay has shown [10] and Jauregui has proved [11] that 
longer codes can be generated from the kernal by certain 
transformations. Given a complementary sequence pair A 
and B with elements a]_ ,a 2 , . . . ,an and b]_ ,b 2 , . . ,bn 

respectively then a longer complementary sequence pair can 
be written as C = AB and D = AB, where AB is a^ , a 2 . . . aj^ ,bi , 
b 2 ,...bj^ and AB is a^^ ,a 2 . . Another opera- 

tion which will produce valid longer complementary sequences 
is the interlacing of the elements. A sequence can be 
written C = A interlace B and D = A interlace B. Interlace 
is defined to be a^^ ,bi ,a 2 ,b 2 , . . . a^^,b^ as required for C and 
D. Once the longer codes have been written by either 
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method then even longer codes can be constructed utili- 
zing the newly created sequence as the given sequence for 
A and B. This operation can be carried on to any desired 
code length. The longer code length will be a function 
of the original code length N. 

It is therefore evident that 2^ (64) valid combinations 
of a kernal can be generated and from any of these comple- 
mentary codes of longer length can be constructed. 

Totally orthogonal complementary*' sequences can be 
generated such that the two pairs of complementary sequences 
when cross-correlated produce only zeros. C. E. French [12] 
in his Master's thesis has shown that every complementary 
code set has exactly one set of two codes that are totally 
orthogonal to it. It is only necessary to know the kernels 
for the original complementary code in order to construct 
the totally orthogonal code. The operation on the comple- 
mentary sequence pair to form the totally orthogonal 
sequence is given by French [12] as follows. Given a 
complementary sequence pair A and B then the totally 
orthogonal sequence pair is B and A. The overbar defines 
the altering of the elements as previously described and 
the underbar defines the inversion of the order of the 
elements. The cross-correlation of A with B added to the 
cross-correlation of B with A will result in the null 
vector (all zeros) . An example of a totally orthogonal 
sequence pair is shown below. 
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Original code .... +1, +1, +1, -1 and +1, +1, -1, +1 

Orthogonal code .. -1, +1, -1, -1 and -1, +1, +1, +1 

B. P. Schweitzer [13] describes generalized complemen- 
tary code sets including non-interacting (totally orthogonal) 
sets in his Doctoral dissertation and extends the code sets 
to include other than binary elements. 

A. G. Fabiila, J. H. Whitehouse, J. M. Speiser and 
others [14,15] have done some preliminary work in an attempt 
to beam form single codes at the Naval Undersea Warfare 
Center, San Diego, California. 

Ultrasonic delay lines have been investigated [16,17, 

18] by J. V. Vollmer, J. H. Whitehouse and J. M. Speiser 
as devices which can be used to physically realize the 
coding of a signal. 
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III. THEORY AND MODELING 



Coding theory has been widely used for the coherent 
transmission and reception of electromagnetic signals to 
enhance certain aspects of communication and radar systems. 
The system generally involves the controlled transmission 
of information. By controlled it is meant that the system 
has control over the transmission as well as the reception 
of the signal. This coherent feature is responsible for 
the success of these coded systems. 

Direction finding is involved only with the reception 
of the signal. The sorting out of the angle-of-arrival 
information has historically been accomplished by uniquely 
defining the amplitude, phase or time for a given AOA. 

These methods were briefly discussed in chapter one. 

A. THEORY OF OPERATION 

The technique presented here involves a combination 
of parameters. The combining, amplifying and altering 
of the individual received signals at each antenna result 
in a unique signal which upon processing provides the 
desired AOA information. 

The type of antenna required is not critical, any 
omnidirectional antenna can be used. The only important 
criterion is that the antennas must be quite similar, if 
not nearly indentical in receiving characteristics. The 
antenna separation should be such that mutual coupling 
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between antennas is held to a minimum. The antennas do 
not form into an array but rather into a corporate set. 

The number of antennas required will be a function of the 
code length. The individual antennas can be any design 
(voltage probes for example) as long as the resultant 
antenna patterns are omnidirectional. A sketch of an 
antenna group for four antennas is shown in figure 3-1. 

The separation between the antennas need not be uniform 
but rather can be any value provided the spacing is known. 
In practice the separation should be as uniform as possible 
to reduce the complexity of the correlation receiver. 

The next functional group after the antennas is the 
signal encoder. The function of this group is to phase 
code each individual antenna output and then to combine 
the individually encoded signals so that this summed 
signal can be applied to a signle transmission line. A 
typical signal encoder for a code of length equal to four 
is depicted in figure 3-2. The signal from each antenna 
is applied to a two output multicoupler type device which 
can provide both a normal (plus in-plus out) output and 
an inverted (plus in-negative out) output. All the normal 
and inverted outputs will be identical except for the delay 
due to antenna separation and the built in 180 degree shift 
separating the normal and the inverted outputs. For a 
signal incident upon the antenna system as shown in figure 
3-1 for an AOA equal to zero degrees (0=zero) the phase 
delay due to antenna separation is equal to the antenna 
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FIGURE 3-1 



30 



SIGNAL ENCODER 



CHANNEL A 




CHANNEL B 



FIGURE 3-2 
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separation. The phase delay due to antenna separation 
will reduce to zero as the AOA of the signal moves off 
and becomes normal to the line of the antennas. The phase 
delay varies from maximum to zero in accordance with the 
cosine of the AOA of the signal. A signal which propagates 
along the line of antennas has an AOA of zero degrees and 
a signal propagating normal to the line of the antennas 
has an AOA of ninety degrees. 

For the example of four antennas there would be eight 
outputs available. The appropriate normal and inverted 
outputs are then combined as required for the complemen- 
tary code in use. For example if the code length is four 
then the required complementary code would be the two 
sequences +1, +1, +1, -1 and +1, +1, -1, +1. The +1 
indicates a normal output and the -1 indicates an inverted 
output. If the antenna nearest the arriving wavefront is 
#1 and the next #2 and so forth, then the signals at the 
antenna side of point A (figure 3-2) would be normal, 
normal, normal, inverted. At the antenna side of point B 
(figure 3-2) the signals would be normal, normal, inverted, 
normal. Since the individual antennas are separated in 
space then the phase separations referenced to antenna #1 
is: antenna #2 one delay, antenna #3 two delays, antenna 

#4 three delays. 

If an incoming signal as shown in figure 3-3 with an 
AOA of zero degrees is assumed then the signals at the 
antenna outputs would be as shown in figure 3-4 and the 
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NOTE; THE UP ARROW ( t ) INDICATES A NORMAL OUTPUT, 

THE DOWN ARROW (*) INDICATES AN INVERTED OUTPUT. 



summed signals at points A and B would be as shown in 
figure 3-5. The coded and summed signals at points A and 
B are now available for transmission from the antenna 
location to a central processing point. A distinct advan- 
tage inherent in this technique is the reduction in the 
number of transmission lines required from the antenna 
location to the processing location. Each of the two 
channels A and B contain all the data required to not 
only determine the AOA but also the information contained 
in the signal. The simultaneous extraction of the DF data 
and the signal itself is somewhat unique. 

The coded and summed signals that is channel A and 
channel B are each applied to banks of tapped delay lines 
in the correlation type receiver. Figure 3-6 depicts a 
typical unit of this section of the system. The terms 
autocorrelation and cross-correlation are used to describe 
the operation depicted in figure 3-6 and delineated in 
the following paragraph. 

The delay lines for the autocorrelation and the cross- 
correlation outputs are identical and there is one set of 
delay lines for each AOA for each channel. The outputs in 
each delay line tap are both a portion of the autocorrela- 
tion output and the cross-correlation output. The 
autocorrelation output is the time inversion of the antenna 
code. For example if the antenna code for a four antenna 
system is +1, +1, +1, -1 and +1, +1, -1, +1 then the 
autocorrelation output is -1, +1, +1, +1 and +1, -1, +1, +1. 
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C/C: CROSS- CORRELATION OUTPUT 
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FIGURE 3-6. 
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The cross-correlation output is the complement of the 
second sequence in the antenna code and a duplication of 
the first sequence in the antenna code. For example for 
the antenna code above the cross-correlation output is 
-1, -1, +1, -1 and +1, +1, +1, -1. The following table 
summarizes the antenna coding and the associated correla- 
tion output. 



Antenna code Autocorrelation Cross-correlation 



Channel A 


+1, 


+1, +1, -1 


-1, +1, +1, 


+1 


"1, “1, +1, 


-1 


Channel B 


+1, 


+1, -1, +1 


+1, -1, +1, 


+1 


+1, +1, +1, 


-1 



The outputs from the autocorrelator and the cross-correlator 
are each summed. The resultant outputs from each correlator 
are shown in figure 3-7 for an exact match. An exact match 
occurs where the delays in the correlators are identical 
to the delays at the individual antennas. The autocorrela- 
tion outputs, A and B are summed directly point by point 
to produce the required signal output. This final output 
as shown in figure 3-8 and termed the autocorrelator 
output exhibits a zero sidelobe level and a mian lobe 
consisting of a duplication of the original input signal 
with an amplitude equal to eight times the original value. 
The multiplying factor of eight is due to the utilization 
of a code of length four since the summed output has a 
value equal to twice the code length. The two cross- 
correlation outputs are also summed directly point by 
point and for this matched condition produce a zero value 
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for all time. This summed signal is shown in figure 3-8 
and is termed the cross-correlator output. 

In summary if the antenna code is composed of elements 
aj^ and (i=l,2,...n) for a code of length n and for j = 
number of possible separations ( j=0 , 1 , 2 , . . . n-1) then the 
autocorrelation series are: 



i=n-j 

Cj= 



i=n-j 

and d ■= I b^bj^^^ 

-i —1 



The final autocorrelation output (figure 3-8) is given by: 

fo j 0 

Cj + dj = J 

1 2n j = 0 



In like manner for a cross-correlation code with the 
above antenna code and j as defined above the final cross- 
correlation output is ej + f j = zero for all values of 
j where ej and fj are the cross-correlation series. The 
cross-correlation series ej and f j are constructed simi- 
larly to the autocorrelation series Cj and d j . 

It can be seen from the preceding that every delay 
in the ninety degrees from an AOA equal to zero degrees 
to an AOA equal to ninety degrees is unique since the 
cosine of the AOA takes on a unique value for all angles 
in this quadrant. 

Signals arriving in the other quadrants 090 to 180, 

180 to 270 and 270 to 000 are processed in a similar manner. 
The antenna code and the correlator outputs are identical 
in the 000 to 090 and the 270 to 000 quadrants therefore 
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the only addition to the system for this region is the 
incorporation of a simple auxiliary antenna group. This 
antenna and associated detectors is used to sort out the 
two quadrants by relative signal strength. The ambiguity 
in this sector is resolved in favor of the stronger signal 
appearing in the auxiliary antenna. The auxiliary antenna 
is required to have two unidirectional patterns physically 
displaced by 180 degrees. Two individual antennas with 
essentially cardioid patterns are adequate for this func- 
tion. The antenna code and the correlator outputs are 
identical for the 090 to 180 and the 180 to 270 quadrants. 
The antenna code utilized in this region is the time 
inverse of the code utilized for the 270-000-090 region 
since the signals arrive at the antenna from the "rear" 
relative to the signals in the 270-000-090 region. An 
additional auxiliary antenna group and detectors are 
required to sort out the two quadrants in this "rear" 
region. 

An ambiguity exists for elevation angles greater than 
zero degrees. The preceding discussion was limited to 
zero elevation angle signals. For non zero elevation 
angle signals the actual delay is the zero elevation angle 
delay multiplied by the cosine of the elevation angle. 

The indicated angle resulting after processing will be in 
error. The error will be such that the value of the 
cosine of the true zero elevation angle AOA would be 
reduced by the cosine of the elevation angle. The resulting 
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angle as processed by the receiver will be increased in 
value from the true AOA for the 000-090 and 180-270 
quadrants and reduced in value for the 090-180 and 270-000 
quadrants. The error is small for low elevation angles 
and also small for angles-of-arrival greater than 10-15 
degrees off the line of the array. 

In operation the outputs of the cross-correlators 
are scanned to find the unique output with the zero value. 
All unmatched cross-correlators will have an output and 
the matched cross-correlator will have no (zero) output. 

The outputs are indexed according to the AOA represented 
by the delay lines therefore the AOA is known as soon as 
the zero cross-correlator output is located. The auto- 
correlator output for the identical delay line is the 
signal of interest which has been increased in signal to 
noise ratio by the correlation gain. Figure 3-9 is a 
simplified block diagram of the complete system. The 
frequency selection circuits of the system included in 
the correlator receiver are the usual devices utilized to 
pass a selected frequency and reject all others. The 
output of the system is therefore unique in two domains. 

The frequency selection section sorts through the frequency 
domain and selects the frequency of interest and the 
correlation processor sorts through the spatial domain 
and a determination of the AOA is made. 

A computer model was developed and utilized to search 
for ambiguities and to determine the limiting factor (s) 
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SIMPLIFIED BLOCK DIAGRAM 



affecting accuracy. The hoped for results of no ambigu- 
ities in a given quadrant was verified and the accuracy 
to which the AOA could be resolved was shown to be 
primarily a function of the delay lines. That is if 
delay lines could be physically realized which were unique 
for each degree of coverage required then the resolution 
of the system is in one degree increments. If delay 
lines stepped in delay by 0.5 degree are realizable then 
the resolution is in 0.5 degree steps and so on for any 
delay line increment. In practice if an AOA fell between 
two delay lines due to the discrete nature of the lines 
then the output to be scanned for would be the minimum. 

In the exact match case this smallest output is equal 
to zero. 

B. SYSTEM MODELING 

The system was modeled utilizing a digital computer 
to provide for flexibility of inputs and also to provide 
a means of obtaining permanent annotated data plots. The 
computer program is discussed in detail in the appendix 
and therefore it will not be detailed in this section. 

Essentially the simulation provides for an input 
signal of any character, any number of antennas, a signal 
encoder, a correlation receiver and a graphical output. 

Any valid complementary code can be utilized. The program 
will generate the proper number of signals and subject to 
the users command will produce either an autocorrelation 
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or a cross-correlation output plotted as amplitude versus 
time. Additionally bounds can be placed on the time range 
over which the correlation takes place. The user can also 
command the program to produce the plotted output over any 
range of angles and stepped in any increment of angle and 
any degree of delay line mismatch. 

Although any code length can be utilized with this 
program a code length equal to sixteen was arrived at as 
a good typical code length. Additionally since the infor- 
mation obtained in one ninety degree quadrant is essentially 
common with the other three ninety degree quadrants the 
investigation was limited to the range of incoming signals 
in the 000-090 degree quadrant. The elevation angle 
utilized was zero degrees. 

1 . Autocorrelation Results 

The incoming signal utilized was a pulse equal 
in time duration to twice the maximum antenna separation 
time delay. The time duration of the correlation process 
was sufficient to encompass all pre-summed sidelobes 
generated by the correlation receiver. The amplitude of 
the incoming signal was taken to be unity. The input 
signal time computation increments were carefully arrived 
at in order to have a priori information to aid in the 
analysis of the data plots. 

The first plot figure 3-10 shows the result of 
autocorrelation with a completely mismatched correlator. 

That is the delay in the correlator is far removed from 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL 



ANGLE OF ARRIVAL* 60.0.. .THIS FILTER IS MATCHED TO 72.6. 



TIME IS PLOTTED ALONG THE X-AXIS ANO THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 ANO 1600. 0/48. OTHEY ARE USEO TO SCALE THE 
PLOT. 



THE FILTER 0ELAY*0.600 TIMES THE ANTENNA SEPARATION 



X-SCALE: "*"* 0.200E 02 UNITS 
Y-SCALE; •'♦"* 0.933E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 
THE WIDTH OF THE MAIN LOBE IS 3.00. 

THE LEAOING EOGE OF THE MAIN LOBE OCCURS AT 371.00. 



FIGURE 3-10 
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the actual matched value. The label above the plot stating 
that the filter delay equals X.XXX times the antenna 
separation indicates the degree of mismatch. The delay 
line filter is matched to the incoming signal when the 
value of X.XXX is 1.000. The angle to which the filter 
(correlator) is matched is also indicated on the printout. 
The next plots figures 3-11 to 3-15 show a much closer 
agreement with the matching of the 'filter to the required 
delay. The expected output for an exact match will be a 
rectangular pulse of original width (duration) and ampli- 
tude equal to twice the code length. In this case the 
width is 100 time units and the amplitude is equal to 32 
amplitude units. Although the output is not "clean" the 
correct output can be seen in the plot. This correct 
output becomes more and more evident as the filter delay 
approaches the matched condition. Figure 3-16 is for 
the matched case. The predicted 100 time units by 32 
amplitude units is the correct output. 

A unique factor concerning this output is that the 
leading edge of the pulse is delayed a specific delay for 
a specific angle-of -arrival . As the angle changes value 
the pulse "walks" towards the left hand ordinate. The 
delay from the left-hand ordinate is equal to the antenna 
delay for the angle-of-arrival times N-1 where N equals 
the number of antennas. Figures 3-17 to 3-21 show this 
shift in the output. Therefore the AOA can be determined 
from the delay of the leading edge of the pulse. While 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL* 60.0. ..THIS FILTER IS MATCHED TO 63.3. 



TIME IS PLOTTEO ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 ANO 1600.0/48. OTHEY ARE USEO TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.900 TIMES THE ANTENNA SEPARATION 



0.0 400.00 

48.00 



X-SCALE; 0.200E 02 UNITS 

Y-SCALE: "♦"* 0.900E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 401.00. 
THE WIDTH OF THE MAIN LOBE IS 68.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 367.00. 



FIGURE 3-11 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0,. .THIS FILTER IS MATCHED TO 62.6. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600. 0/^8. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.920 TIMES THE ANTENNA SEPARATION 



« 4t4< 4i * * * ** *4- « 






X-SCALE: "*"= 0. 200E 02 UNITS 
Y-SCALE: "♦••= 0.933E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOSE IS A08.00. 
THE WIDTH OF THE MAIN LOSE IS 74.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 371.00. 



FIGURE 3-12 
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the following data is for a new angle of arrival. 



ANGLE OF ARRIVAL* 60.0. ..THIS FILTER IS MATCHED TO 62.0. 



Pl6t^°’° ° 1600.0M8.0THEY ARE USED TO SCALE THE 



46.00 



THE FILTER DELAY=0.940 TIMES THE ANTENNA SEPARATION 



400.00 800.00 1200.00 1600.00 

X-SCALE: "♦"= 0.200E 02 UNITS 
Y-SCALE: ••*••= 0.933E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 408.00. 
THE WIDTH OF THE MAIN LOBE IS 74.00. 

THE LEADING EDGE OF THE MAIN LG6E OCCURS AT 371.00. 



FIGURE 3-13 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0, ..THIS FILTER IS MATCHED TO 61.3. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600. 0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY-0.960 TIMES THE ANTENNA SEPARATION 




■ 6.00 . + 
**-t-*1ftL*^f****^iii**1tf**ti**-t-»********-t-****:ti****ii-*i********^*ti****ti**i-***ti*****-^*********H-** 

0.0 400.00 800.00 1200.00 1600.00 



X-SCALE: 0.200E 02 UNITS 

Y-SCALES 0.900E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 417.50. 
THE WIDTH OF THE MAIN LOBE IS 83.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 376.00. 



FIGURE 3-14 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL^ 60.0. ..THIS FILTER IS MATCHED TO 60.7. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.980 TIMES THE ANTENNA SEPARATION 






12.00 



3.00 




1filfi.tlH^^t:titf*t^ilf4.^^^****t^ilf^t‘****************************^*********-t-**ti******-^*********+** 

0.0 400.00 800.00 1200.00 1600.00 



X-SCALE: 0.200E 02 UNITS 

Y- SCALE; 0.900E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 417.50. 
THE WIDTH OF THE MAIN LOBE IS 83.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 376.00. 



FIGURE 3-15 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL* 60.0. ..THIS FILTER IS HATCHED TO 60.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600. 0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=1.000 TIMES THE ANTENNA SEPARATION 



0.0 400.00 800.00 1200.00 1600.00 

m ***** *-t- rn****+tf***** ***+****** ***^*********^*******ti#.^.*^**i^m***^ ****** ***■*■** 

48.00 ♦ . + 

* * 



0.0 




************************************************************ **-i. ************ *******+** 

0.0 400.00 800.00 1200.00 1600.00 



X-SCALE: 0.200E 02 UNITS 

Y-SCALE: 0.800E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 425.00. 
THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 375.00. 



FIGURE 3-16 



THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 50.0. ..THIS FILTER IS MATCHED TO 50.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=1.000 TIMES THE ANTENNA SEPARATION 



X-SCALE: ••♦«= 0.200E 02 UNITS 
Y-SCALE; «♦"= 0.800E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 
THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 479.99. 



FIGURE 3-17 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 



ANGLE OF ARRIVAL* 60.0. ..THIS FILTER IS MATCHED TO 60.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

°0. 0/0.0 AND 1600.0/<>8. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER OELAY=1.000 TIMES THE ANTENNA SEPARATION 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*******+******°**°*******************°*********** 

A8.00 ♦ 



iiii*i+****i*iii4*******************+'^************** 
^*00. 00 800.00 1200.00 



X-SCALE: 0.200E 02 UNITS 

Y-SCALE: ”♦”= 0. 800E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS ^25.00. 
THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 375.00. 



FIGURE 3-18 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL*: 70.0. ..THIS FILTER IS MATCHED TO 70.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=1.000 TIMES THE ANTENNA SEPARATION 

0.0 400.00 800.00 1200.00 

48.00 * 



8.00 



0.0 




*41 4^41* 4t4i *■(-«** *4. 4t*<i*««***44t«**«4t***4#4t**;4*««*4^* 

0.0 400.00 800.00 1200.00 1600.00 



X-SCALE: "♦"= 0.200E 02 UNITS 
Y-SCALE; "*«= 0.800E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 305.00. 
THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 255.00. 



FIGURE 3-19 
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THE FOLLOWING OATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 80.0. ..THIS FILTER IS MATCHED TO 80.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND l600.0/<i8. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER OELAY=l.OOO TIMES THE ANTENNA SEPARATION 



X-jSCALE: 0.200E 02 UNITS 

Y-SCALE: 0. 800E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 170.00. 
THE WIDTH OF THE MAIN LOBE IS 100.00, 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 120.00. 



FIGURE 3-20 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL* 90.0. ..THIS FILTER IS MATCHED TO 90.0. 



THE FILTER 0ELAY=1.000 TIMES THE ANTENNA SEPARATION 



0.0 + * 

0.0 400.00 800.00 1200.00 1600.00 

X-SCALE: ••*’•* 0.200E 02 UNITS 
Y-SCALE: ••*••= 0.800E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 50.00. 

THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-21 
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this is a true statement it is not the criterion utilized 



for AOA determination. The disadvantage of this method 
of establishing the AOA is the requirement that the signal 
have a sharply discernible modulation event. An unmatched 
autocorrelation output will usually have an output which 
exceeds twice the code length, 32 in this case, at some 
time in the correlation process. This was also rejected 
as a method of determining the exact match case since this 
event generally had an extremely short time duration. 
Therefore the autocorrelation output is only utilized to 
provide information concerning the signal and the associ- 
ated modulation. The cross-correlation output establishes 
the AOA and is required to indicate which autocorrelation 
output is correct for the true determined AOA. 

2 . Cross-Correlation Results 

The program was processed utilizing the signal and 
the system parameters discussed for the autocorrelation 
process. The changing of one control card in the data 
deck is all that is required to shift the program from 
autocorrelation to cross-correlation. Figure 3-22 is a 
plot of a completely mismatched case. The expected output 
for the matched case is zero for all time. As can be seen 
from the plot for this case the output is not zero. Figure 
3-23 to 3-27 are output from conditions where the delay is 
more closely aligned to the required matched case. Even 
though the filter (correlator) is more closely matched to 
the correct delay a null value cannot be seen to be rapidly 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0. ..THIS FILTER IS MATCHED TO 72.6. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS.. DISREGARD 

0. 0/0.0 AND 1600.0/A8. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.600 TIMES THE ANTENNA SEPARATION 



X-SCALE: "*"= 0.200E 02 UNITS 
Y-SCALE: «*•«= 0.127E 01 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 0.0 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-22 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0.. .THIS FILTER IS MATCHED TO 63.3. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/<»8.0THEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAYED. 900 TIMES THE ANTENNA SEPARATION 



-lA.OO ♦ 

0.0 



X-SCALE: 

Y-SCALE: 



0.200E 02 UNITS 
•»♦«= 0.103E 01 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-23 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0. ..THIS FILTER IS MATCHEO TO 62.6. 



TIME IS PLOTTED ALONG THE X-AXIS ANO THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/A8.0THEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.920 TIMES THE ANTENNA SEPARATION 



1600.00 

. + 



- 12.00 + 

0.0 AOO.OO 800.00 

X-SCALE: '•♦'•= 0.200E 02 UNITS 
Y-SCALE: "♦"= O.IOOE 01 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0. 



FIGURE 3->24 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0. ..THIS FILTER IS MATCHED TO 62.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600. 0/^8. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER 0ELAY=0.940 TIMES THE ANTENNA SEPARATION 






X-SCALE: ”♦"= 0.200E 02 UNITS 

Y-SCALE: O.IOOE 01 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-25 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0.. .THIS FILTER IS MATCHED TO 61.3. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.960 TIMES THE ANTENNA SEPARATION 






X-SCALE: 0. 200E 02 UNITS 

Y-SCALE: O.IOOE 01 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 0.0 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-26 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0.. .THIS FILTER IS MATCHED TO 60.7. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 
OUTPUT IS PLOTTED ALONG THE Y-AXI S. ...... .DI SREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=0.980 TIMES THE ANTENNA SEPARATION 



0.0 

48.00 + 




12.00 + . ♦ 



0.0 


400.00 




800.00 


1200.00 


1600.00 


X-SCALE: 

Y-SCALE: 


”♦"= 0.200E 02 UNITS 
0. lOOE 01 UNITS 












THE TIME OF OCCURANCE OF THE CENTER OF THE 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 


MAIN LOBE IS 

0.0 . 


0.0 . 



FIGURE 3-27 
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approaching as the delay approaches the correct value. 

The correct value is indicated on the plot for X.XXX 
equal to 1.000 in the statement just above the plot. The 
exact case is shown in figure 3-28. This zero output is 
so unique and dramatic that it was chosen as the criterion 
to establish the exact AOA. 

Since the autocorrelation and the cross-correlation 
outputs are processed from the same set of tapped delay 
lines for each AOA then as soon as the zero output is 
located the autocorrelation output of the signal is 
available for further processing and/or analysis. 



67 



THE FOLLOWING DATA IS FOK A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL* 60.0. ..THIS FILTER IS HATCHED TO 60.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY*1.000 TIMES THE ANTENNA SEPARATION 



^*^^^4.^*^i^^ti*t‘*-t-*^***m****-*********4-************ 
800.00 1200.00 1600.00 



X-SCALE: 0.200E 02 UNITS 

Y-SCALEs 0. 800E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 0.0 
THE WIDTH OF THE MAIN LOBE IS 0.0 . 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 0.0 . 



FIGURE 3-28 
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IV. SUMMARY AND RECOMMENDATIONS 



A. SUMMARY 

The problem of determining the angle-of-arrival of a 
propagating radio wave has been investigated from the 
viewpoint of simplifying the antenna system required 
while retaining the inherent accuracy of a complex antenna 
system. The technique arrived at utilized the phase 
coding of the received signal and the processing of this 
signal for two outputs. The required outputs are a pseudo 
autocorrelation output and a cross-correlation output. 

The realization of AOA information and the simultaneous 
ability to recover an autocorrelation type output of the 
signal is detailed in chapter III. In addition to the DF 
and the intercept information this technique also results 
in a simple configuration from the antenna location to the 
processor, namely two transmission lines vice the multiple 
lines required by other systems. 

The accuracy of a DF system has historically resided 
in the ability to produce either a large complex antenna 
system or the wide separation of simple antennas or in 
some cases a little of each. The technique described in 
this investigation requires simple antennas in a simple 
configuration. Additionally the antenna system is 
stationary . 
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The number of antennas is determined by the code length 
desired. The code length is chosen in accordance with the 
desired correlation gain and the physical environment. 

The correlation gain is equal to twice the code length. 

The bearing accuracy or resolution ability is a direct 
function of the incremental steps in the delay lines. 

That is if the delay lines are matched in delay to one 
degree increments then the resolution of the system is 
one degree. If the delay lines are in increments of half 
a degree then the system will resolve lines of bearing to 
half a degree. Therefore the resolution ability is limited 
by the ability to physically produce these lines in the 
increments desired. 

This technique will produce ambiguous AOA data since 
the system can only uniquely determine AOA data in a 180 
degree sector. The data in the mirror image 180 degree 
sector can be confirmed (or discarded) by a simple 
auxiliary antenna sampling signal strength. Therefore 
the ambiguity can be resolved. 

A digital computer model was utilized to show that 
there are no ambiguous angle-of-arrivals in the 180 degree 
sector. The model also verified that the resolving ability 
of the system, was only a function of the ability to 
physically realize the appropriate delay lines. 

B. CONCLUSIONS 

The ability to determine the angle-of~arrival of an 
incoming signal was shown to be a processing function 
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utilizing a relatively simple antenna arrangement. 
Ambiguities could be resolved and the resolution ability 
of the system was not limited by the theory presented but 
rather is a function of physically realizable hardware. 

C. RECOMMENDATIONS FOR FURTHER INVESTIGATION 

1. A simple antenna system, signal encoder, delay 
lines and receiver should be built and data taken to verify 
the computer model results. A four antenna group, simple 
encoder, delay lines for thirty degree steps and any 
available receiver could be utilized. The frequency 
chosen should be low enough such that phase lengths are 
not a problem, but not so low as to require extreme 
spacing between antennas in order to reduce mutual coupling. 
A frequency in the region of 200 to 300 MHz would probably 
suffice . 

2. A study should be made to attempt to reduce the 
number and/or precision required for the tapped delay 
lines . 

3. A study should be made to determine parameters 
required for the optimum performance utilizing the simplest 
antenna. The same study should undertake the determination 
of the optimum performance to be obtained by utilizing more 
than one simple antenna system in various configurations. 

4. An investigation should be made of a means to 
negate the error due to a non-zero elevation angle of an 
incoming signal. 
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APPENDIX 



COMPUTER PROGRAM 



The computer simulation is a digital program written 
in FORTRAN IV level G. It was processed by an IBM 360/67 
system and the built in functions are those normally 
available on this system. 

The program requires a data deck input for the specific 
complementary code used in addition to the parameters 
required to process the code. The following words are 
peculiar to this program. 

CODEl: First sequence in a complementary code. 

CODE2: Second sequence in a complementary code. 

SIGN: Summation sign, used for addition of correlator 

outputs 

DELDEG: Incremental angle shift in degrees, used to 

obtain a new value for the angle. 

START: The starting value in degrees of the incoming 

signal . 

STOP: The final value in degrees of the incoming 

signal . 

START-T: The starting time for the computation. 

STOP-T: The completion time for the computation. 

DSTEP: Incremental step for the parameter D, used to 

obtain a new value for the filter (correlator) 
delay. 
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NRANTS : 


The number of antennas in use. 


SPACE : 


Spacing between antennas; antennas have 
uniform spacing. 


DELT: 


Computation time increment. 


TIMES: 


Multiplying factor used in the subprograms 
in order to utilize integer arithmetic and 
logic . 


JJ: 


Control for surpressing data printout; do not 
print if JJ=2. 


JJJ: 


Control for correlation function; if JJJ=1 
then cross-correlate, if not then autocorrelate 


JJJJ: 


Control to surpress the final plot which is a 
summary of maximum amplitude versus angle; if 
JJJJ=2 then do not plot. 



The data deck must be assembled in a specific sequence 
in order to have the program operate properly. Addition- 
ally if the data is not known a card must still be loaded 
into the data deck in order not to invalidate the read 
instructions. The data must be provided whether valid 
or not. The data cards are assembled as in figure A-1. 

The first card in the data deck is CODEl and the last is 
JJJJ. All data cards are punched starting in column one. 
CODEl and C0DE2 are formated floating point F4 . 1 and each 
card can have up to eighteen entries. Each code is read 
from one card, if longer codes (greater than 18) are 
required then the read statement must be changed to read 
more than one card for each code. NRANTS, JJ, JJJ and 
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TYPICAL DATA DECK 
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FIGURE A-1 
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JJJJ are formated in integer 111 all other data cards are 
formated in floating point F11.6. 

The program is set up to provide for a pulse signal 
input with the pulse duration equal to twice the delay 
due to antenna separation. The program will automatically 
scale the built in pulse duration (W) and the antenna 
separation (SPACE) to accommodate any code length. 
Additionally the computation time is automatically scaled 
in order to include all the correlation outputs greater 
than zero and also provide for the maximum number of 
computations per angle and filter delay setting. The 
maximum number of computation intervals for a fixed angle 
and delay is 2000. The codes are read into computer 
memory and all necessary operations are accomplished in 
order to provide the proper autocorrelation and cross- 
correlation functions as a function of time, angle and 
filter delay. If data is provided which is beyond the 
array storage set aside for the computation then the 
program will attempt adjustments in the size of the data. 
After twenty attempts with no success the program will 
terminate and print out: DATA EXCEEDS ARRAY STORAGE... 

JOB TERMINATED. 

The program will input all parameters, compute other 
required parameters in accordance with the inputted param- 
eters, compute an output for a fixed time, angle and delay 
and then step the time variable to the next increment. 
After the time variable has been stepped across its range 
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the angle variable is incremented to a new value and the 
output is computed for all time values as the time vari- 
able is again swept through its range. After the angle 
variable has been swept through its range the delay 
variable is incremented to a new value. The output is 
then computed across the time range and the angle range. 
This sequence is followed until the full range of values 
for the delay variable have been used. In this manner 
the outputs of the filters (correlator outputs) are 
computed for all times of interest, all angles of interest 
and all desired variations of the filter delay lines. 

The output available from this program is a plot 
routine of amplitude versus time for each and every value 
of angle and filter delay required. Additionally as shown 
in figure A-2 and A-3 several important statements are 
printed. The statements below the plot are valid and 
useful only for the matched autocorrelation condition 
where D=1.00. The data in these statements can be used 
for the autocorrelation unmatched case to roughly locate 
the pulse shape as the variable D approaches unity. After 
all values have been computed for the full range of time 
and angle variables and prior to the stepping of the 
variable D a summary plot of maximum amplitude versus 
angle and a table of this data can be plotted by selection 
of the proper control card. 

A listing of important parameters and statements 
concerning the interpretation of the various control 
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commands are printed at the beginning of each computer 
run to aid in analysis of the data and the plots. 

The maximum core size is 150K and due to the fine 
grain computation required by the program execution time 
for most investigations is in the 15 to 20 minute region. 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL= 60.0. ..THIS FILTER IS MATCHED TO 63.3. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/48. OTHEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAYED. 900 TIMES THE ANTENNA SEPARATION 



4- 4c # 4t 4t « « « « 4= + « « « « * 



» 4c 4c # > 4c 4c 4c 4c 4c A * 4c 4c 4c * 4c 4c 4c 4c 4c * 4c -f « « * * * * * * -I- 4' 41 * -i- ♦ # 
800.00 1200.00 1600.00 



X-SCALE: 
Y- SC ale: 



"4c"= 0.200E 02 UNITS 
"♦"= 0.900E 00 UNITS 



THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 401.00. 
THE WIDTH OF THE MAIN LOBE IS 68.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 367.00. 



FIGURE A-2 
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THE FOLLOWING DATA IS FOR A NEW ANGLE OF ARRIVAL. 

ANGLE OF ARRIVAL* 60.0.. .THIS FILTER IS MATCHEO TO 60.0. 



TIME IS PLOTTED ALONG THE X-AXIS AND THE FILTER 

OUTPUT IS PLOTTED ALONG THE Y-AXIS DISREGARD 

0. 0/0.0 AND 1600.0/A8.0THEY ARE USED TO SCALE THE 
PLOT. 



THE FILTER DELAY=1.000 TIMES THE ANTENNA SEPARATION 

0.0 400.00 800.00 1200.00 

48.00 + 



0.0 + + 

♦ ♦ + ♦*«<=*♦♦** + ****♦♦#♦♦ + **♦* ♦♦♦♦♦+«***.*4:** 

0.0 400.00 800.00 1200.00 1600.00 

X-SCALE: "*"= 0.200E 02 UNITS 
Y-SCALE: "*"= 0.800E 00 UNITS 

THE TIME OF OCCURANCE OF THE CENTER OF THE MAIN LOBE IS 425.00. 

THE WIDTH OF THE MAIN LOBE IS 100.00. 

THE LEADING EDGE OF THE MAIN LOBE OCCURS AT 375.00. 



FIGURE A-3 
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PHASE CODED DIRECTION FINDING SYSTEM SIMULATION UTILIZING COMPLEMENTARY CODES 
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100 FORMAT (F11.6) 

READ (5,105) NRANTS,JJ,JJJ,JJJJ 
105 FORMAT (111) 
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DELT, TIMES, SI GN, OELOEE,W,DtD STEP, STARR, STOO, SPACE, 
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START COMPUTATION WITH INITIAL OR NEW VALUE OF ANGLE OF ARRIVAL 
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COMPUTE SIGNAL AT OUTPUT OF EACH INDIVIDUAL ANTENNA 
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REVERSE THE CODE 
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